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1 Introduction 



The one dimensional periodic potentials play an essential role in physics. The simple 
model of Kronig-Penney has been recently used for describing superlattice structure of 
ultrathin epitaxial layers with diverse constituent-semiconductor compositions such as 
GaAs and Al z Gai_ x As pQ. The Kronig-Penney lattice with R matrix interaction at the 
end of each block has been proposed for analyzing the behavior of neutrons in crystals 2 . 
The Mathieu equation has successfully described the diffraction of intense standing light 
by strong sinusoidal media [H]. 

The progress of these models is delayed by the fact that few analytically solvable 
periodic potentials are known, including the text-book examples of Kronig-Penney [I], 
Lame, and Mathieu [Ej. The inverse Sturm-Liouville methods could be applied [Hj but 
up to the authors knowledge no simple periodic potentials have been found in this way. 
In order to enlarge the class, the Darboux transformations of various orders can be used 
[7llH|ini[TTHlllj. the subject commonly known as the super •symmetric quantum mechanics 
(SUSY QM) jl2H3llTl] . By employing 2-nd order irreducible techniques, proposed initially 
for non-periodic potentials JE] (see also a hint by Krein JE]), the factorization energies 
above the bound states can give new non-singular potentials |17U18llT§] . We shall now 
extend the method by showing that it can be applied in a non-singular way to the periodic 
potentials if the -E-values are in the same spectral gap. In the most interesting cases, the 
operation can generate the defects of the periodic structure, creating the bound states 
inserted into spectral gaps (a typical phenomenon in solid state physics 20 ). 

Below, the technique will be applied to Lame and periodically continued soliton po- 
tentials. In its most inconspicuous form, it produces an interesting non-local effect: the 
transformed potential becomes an exact or approximate displaced copy of the initial one, 
a phenomenon which we call the translational invariance with respect to Darboux trans- 
formations or shortly Darboux invariance |21j . Quite significantly, the effect shows itself 
asymptotically even if the periodicity of the initial potential is not preserved. It permits 
to understand the structure of Darboux generated lattice impurities as the contact effects 
caused by a conflict between two non-local S USY transformations. 



2 First and second order Darboux 
transformations 

Let us outline briefly the main points of the Darboux method. Consider the Schrodinger 
equation with an arbitrary potential Vq(x): 

h ip E (x) = Etp E (x), h = -d 2 x + V (x). (1) 

The method permits to use the solutions iPe( x ) of the initial Schrodinger equation (^Q) to 
obtain the solutions (Pe(%) of the transformed Schrodinger equation 

h 1 i PE = EipE, h 1 = -d 2 x + V 1 {x) (2) 

by applying a certain differential operator L, (Pe( x ) = LiPe(x)- In the first order case 
L = —d x + w(x), the transformation yields: 

AV(x) = Vi{x) - V {x) = -2 [\nu(x)]" (3) 
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where u(x) is a solution of Q called the transformation function while w(x) = [lnit(x)]' 
is the superpotential. Note that the method applies for any E € M in the spectrum or 
in the resolvent set of ho; it links the solutions of (Q) and ©, without requiring that ipg 
should be bounded or square integrable (though in some problems, such assumptions can 
be pertinent). Quite obviously, in order to avoid the creation of singularities in w(x) and 
in V\{x) one must look for nodeless u(x). 

In the second order case, the new potentials and eigenfunctions are defined by a pair 
of transformation functions u\{x) and u 2 (x), /io^i,2 = 01,2^1,2 (we use deliberately the 
symbols a±, a 2 instead of E±, E2 to stress that both parameters do not need to belong 
to the spectrum of ho). The transformation reads: 

(p E = W~ 1 (u 1 ,U 2 )W(ui,U 2 ,'tpE), f.s 

AV= -2[]nW( Ul ,u 2 )}" { ' 

where W denotes the Wronskian of the corresponding functions. 

For the periodic potentials Vq{x) = Vq(x + T) the method is often affected by singu- 
larities. Some authors opt to use the band edge solutions [71I51IT7] (which are periodic or 
antiperiodic [22] ). We are going to show that the use of the generalized Block functions 
belonging to the same spectral gap brings even more interesting results. 



Bloch eigenfunctions and Darboux 
transformations 



The Bloch functions are usually defined as physically interpretable eigenfunctions of (0) 
for E belonging to the spectrum of ho- Yet, they exist also out of the spectral area. By 
writing for any E £ R as the first order differential equation for a vector formed by tp 
and its derivative ip' 

(5) 

one sees that tp(x), ij)'{x) are given by a certain linear transformation applied to tp(0), ip'(0): 
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where b(x) is a 2 x 2 simplectic transfer matrix j22U23j . Note that b(x) can be expressed 
in terms of special real solutions vi(x), V2(x) of © defined by the initial data vi(0) = 1, 
v' 1 (0)=0, v 2 (0)=0, v' 2 (0) = l: 



b(x) 



vi(x) v 2 (x) 
v[(x) v' 2 (x) 



(7) 



implying Det[6(x)] = 1. If Vq(x) is periodic, an essential role belongs to the Floquet (or 
monodromy) matrix b(T) (see e.g. |22U24| ). Since Det[6(T)] = 1, its eigenvalues f3 are 
given by: 

(3 2 — Dp + 1 = (8) 



where 



D = D(E) = Tr[6(T)] 



0) 
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is called the Lyapunov function, Hill determinant or discriminant of (0) 6l l22ll2"o] . As 
follows from ©, /? takes two values (3± such that /?+/?_ = 1: 

P± = D/2 ± V£> 2 /4 - 1 • (10) 

Whenever for any E£l, one of the eigenvectors of b(T) (for either (3 = /3+ or f3 = /?_) 
is used as an initial condition for ip and ip' at x = 0, it originates a special solution of Q 
for which the transfer law © at x = T reduces to: 

V>(T)=/?V(0), V'(T) = W'(0) (11) 

and more generally 

ip(x + nT) = (3 n i((x), i//(x + nT) = (3 n ip'(x), (12) 

n = 0, ±1, . . . The eigenfunctions ip{x) of Q which fulfill (|11I12|) exist for any parameter 
E = a € R (not necessarily belonging to the energy spectrum of ho) and are called the 
Block functions. Their structure depends essentially on the values of E. 

If \D(E)\ < 2, E is inside of a spectral band and the eigenvalues 0± in (fTOjl are complex 
numbers of modulus 1. We can put then 

(3 + = exp(ikT) /?_ = exp(-ikT) (13) 

where A; is a real parameter called the crystal quasimomentum. The equations (|10j) and 
(|13|) define an implicit function k = k(E) called the dispersion law. The corresponding 
bounded and essentially complex Bloch functions are the ones traditionally considered in 
the solid state physics. 

The equation \D(E)\ = 2, in turn, defines the band edges [22] • Following |T2j let us 
denote them by 

E < E x < E v < E 2 < E 2 > < . . . < Ej <E f <... 

At each band edge the Bloch eigenvalues (fTflj) are /3+ = (3- = ±1, and the degenerate b(T) 
defines one Bloch function, periodic or antiperiodic. Both edge eigenfunctions rpj and ijjji 
are real and have the same number j of nodes. 

If \D(E)\ > 2, the eigenvalues {TDJ) are real (3 + = f3, (5- = 1/(3 (0 ^ G R); b(T) has 
the real eigenvectors for both (3 = /3±, originating the real Bloch functions, which form 
a natural basis for the general solutions of (^Q). (In what follows, when speaking about 
the Bloch functions in this regime, we shall always have in mind the real Bloch functions, 
without specially mentioning it. The complex Bloch functions |26j have a constant phase 
and are easily reduced to the real ones.) Note that all solutions of Q diverge in either 
+oo or -co, hence E is in a forbidden band (resolvent set of ho). The Bloch functions here 
are deprived of a physical meaning but happen to be crucial as transformation functions 
in the Darboux algorithms. Their properties still depend on the specific localization of E 
in the resolvent set. 

The interval (— oo, Eo) constitutes the lowest forbidden band where D > 2 and (3± > 0. 
In the following spectral gaps (Ej,Ejr) one has either D > 2 and (3± > for j even or 
D < -2 and j3± < for j odd. 
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Let u(x) be one of the Bloch eigenfunctions for E = a < Eq. Due to (|12|). u(x) has 
either no zeroes at all or a periodically distributed infinite set. The last possibility cannot 
occur since then, according to the Sturm oscillator theorem (see e.g. |27|). the ground 
state eigenfunction should have an infinity of zeroes as well, which contradicts the fact 
that i/jq(x) is nodeless. Thus, for any a < Eq there exist two linearly independent nodeless 
Bloch eigenfunctions, both suitable to generate non-singular Darboux transformations 
As nodeless must be the Wronskian W (111,112) for any two Bloch functions u\, ui with 
ol\,oli < Eq, thus generating a non-singular reducible transformation Q |28j . 

Let us now examine the situation in the higher spectral gaps. We shall show that while 
the first order Darboux transformations here are singular, the second order ones can be 
regular, in analogy with j^H]. Indeed one has: 

Proposition 1. Let u±, U2 be two linearly independent Bloch eigenfunctions of (fTjl. 
Ui(x+T) = (3iUi(x), JiQUi = ctjiij, with a%, «2 in the spectral gap (Ej, Ej/) of a continuous, 
periodic Vq(x). Then: (i) In each periodicity interval [xo,xo + T) both u\(x) and U2(x) 
have j roots, depending continuously on ct% and Q2; (h) if ai 7^ 0,2, or a\ = 012 = ol but 
01 7^ 02 y then the roots of u\ and U2 cannot coincide; they form two infinite alternating 
sequences extending from —00 to +00; (iii) the Wronskian W(x) = W(u±, U2) has no nodal 
points on R (and so, defines a non-singular transformation 

Proof, (i) is well known in the theory of Hill's equations [221. (h) If ct\ < 0:2 then due 
to the oscillatory theorem [27j, between each two neighbouring nodes of u\ there is at least 
one nodal point of U2 (but no more than one because of (i)). Moreover, the nodal point of 
u\ cannot be nodal for 112, since then the oscillatory theorem would require still one more 
node of 112 in between the vicinal nodes of u±; so 112 would have more nodes than u\ in 
[xq,xq + T). If a\ = o<2 but 0\ 7^ 02, u%, 112 cannot have a common nodal point. If they 
did, they would satisfy proportional initial conditions at the common node and so, they 
would be proportional all over R, which is impossible if 0\ ^ 02- Finally, due to (|12|) the 
zeroes of u\ and U2 are periodically distributed on R, so they must form two alternating 
sequences extending from —00 to +00. (iii) Let now v and v be two neighbouring nodes 
of u\(x), U2(x) respectively. Observe, that at the extremes of the interval (y, v) the values 
W(v) = —U2{v)u' 1 (y) and W(u) = u\(D)u' 2 (v) cannot vanish and must have the same sign. 
Indeed, since we assume that neither u\ nor v>2 have zeros inside of (y, v) the sign of u\(x) 
coincides with the sign of u'Av) while the sign of 1*2 (x) is opposite to that of u' 2 (v) in all 
\u, v\. Moreover, the first derivative W'(u\, U2) = (et\ — ct2)u\U2 either vanishes everywhere 
(if «i = 02), or nowhere in (v,D) (if a\ ^ 02)- In both cases W(u\,v,2) is monotonic, 
extends between two non-zero values of the same sign and so, it cannot vanish in [u, v\. 
Since the same holds for any other neighbouring nodes of u\ and U2 then W(u\,ui) has 
no nodal points on the entire R. I 

While for a fixed a £ (Ej,Eji) the zeroes of the Bloch pair are isolated points, when 
a varies, they draw a sequence of continuous nodal curves. These curves cannot intersect 
inside of (Ej,Eji) and, moreover, each one can intersect any vertical line x = c only 
once (compare Proposition Iii), so they form a sequence of strictly monotonic branches 
which meet at the band edges, where two Bloch solutions degenerate to one (see Fig. 1); 
a pattern which permits to extend the Proposition 1 to non-trivial linear combinations of 
the Bloch functions. In what follows we shall use the simplified symbol vP to denote the 
Bloch eigenfunction for any a, corresponding to the Bloch eigenvalue (so to be exact 
u@ = u(x, a, 0)). 
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Figure 1: The "nodal curves" illustrating the a-dependence of the nodal points of two Bloch 
functions u 13 (black line) and u 1 ^ (shadow line). Notice that any two vicinal nodes situated at 
a = ai are never separated just by a single node at a — a.2, thought they can be separated by a 
pair. The graphic was obtained for the Lame potential (|23[) with n = 1, m = 0.99. 

Theorem 1. Let u\{x) = u l3l (x), U2(x) = u^ 2 {x) be two nontrivial Bloch functions of 
ho for two different eigenvalues ai,«2 € (Ej,Ej/) and let v\(x) = u^ 1 + niu 1 ^ 1 , V2(x) = 
u^ 2 + K2U 1 ^ 2 , where (^1,^2) E M 2 . Then there exists an infinite closed sector in R 2 , 
bordered by n\- and K2-axes, where the Wronskians W(v\,V2) are nodeless. 

Proof. Note that the Bloch functions u^ 1 ,?/ 1 /^ 1 (a = a\) and u^^u 1 ^ 2 (a = 02) 
determine four sequences of nodal points on K. In view of the nodal configuration of Fig. 1, 
each two neighbouring nodes v\,v\ of u^ 1 and u 1 ^ 1 can be separated by a pair of nodes 
U2-, V2 of u@ 2 , u 1 /^ 2 but not by just one of them. Without loosing generality one can choose 
the point x = so that vP' (0) ^ 0, -u 1 /^ (0) ^ and that the four subsequent nodal points 
are v\ , v\ , v 2 , i>i (belonging to u' 31 , u 1 ^ 1 , u 132 , u 1 ^ 2 respectively). By multiplying each 
Bloch functions by ±1 one can also achieve u^{xq) > and u 1// ^ ,; (xo) > 0. According to 
Proposition 1, all Wronskians W(u^,u^), W(uWi, W(itVft jU VA) 

have constant signs; moreover, it is straightforward to check that they are now strictly 
positive. Henceforth, by choosing ki,K2 > one obtains 

W{v u v 2 ) = W(u pl ,u^) + k 2 W(u i3i ,u 1/132 ) + k x W{u x ' p \u^) + KiK 2 VF(n 1 / /3l ,u 1 / /32 ) 

strictly positive, producing a non-singular Darboux transformation ©. I 



4 The translational effects of the Darboux 
operations 

We shall show now that for certain classes of potentials the Darboux transformations 
might yield an interesting non-local effect. Consider first of all the 1-soliton well |3Uj 

V Q {x) = -2 7 2 sech 2 7 x , 70 > 0, (14) 

with one eigenvalue Eq = —79- We shall check that the wells (|14j) admit special Darboux 
transformations leading to the exact coordinate displacements. Note that the Schrodinger 
equation Q with the potential (|14|) can be exactly solved for any E £ M. by applying the 
Darboux transformation (J3J) to the free Hamiltonian. One of the solutions is: 

= cosh7o(x + ^) e _ 7lX 
cosh 7 ox 



6 



If now ai = — 7 2 < Eq, then using (|15|) as a 1-SUSY transformation function for our 
1-soliton well Vq(x), one obtains a displaced version of (fTl]) . Ui(x) = Vo(x + <5i), with 

^! = — artanh— (16) 
70 7i 

If ot\ > Eq, the first order Darboux transformation can neither produce the displaced, 
nor even nonsingular potential, but the second order SUSY opens new possibilities. Taking 
5i = 5 + = 5 + ir', where t' is half the imaginary period of the transparent well (|14j) . 
we easily induce a complex displacement generated by u\ = e~ llX sinh7o(x + 5)/ cosh 702; 
and leading to the real but singular V\{x) = 27qCsc1i 2 7ox. By repeating now the operation 
with a new complex 82 = 5' — i-^- one returns to the original transparent well displaced 
additionally by 5" = 5± + 82 = 8 + 8' . The corresponding nodeless Wronskian in @ is: 



„-, . / cosh7o(x + (5") _ (7l+72 w 

W{ui,u 2 ) = (71 -72) r e w^ 2 > . (17) 

cosh 70 x 

Until now our second order displacements (jlj) were backed by pairs of first order trans- 
formations © of the displacement type, with real or complex <5's. It would be interesting 
to achieve the same effect without this kind of first order scenario. This indeed happens 
for the symmetric 2-soliton well 



= 2(7? ~ 7l)(7isech 2 7ix + 7 2 csch 2 7 2 :r) 
(71 tanh jix — 72 coth 72a;) 2 



(18) 



obtained from the null potential using @ with u\(x) = COSI171X, u 2 (x) = sinh72X, where 
72 > 7i- The potential (fl8|) has two discrete energy levels at Eq = — 7 2 , E\ = — 7 2 . 
Denote now u^(x) = e~ 13X : U4,(x) = e~ 1AX and Wij(x) = W (ui(x) , Uj (x)) , Wiji(x) = 
W(ui(x),Uj(x),ui(x)); then apply an isospectral second order Darboux transformation 
to the potential (|18jl using two new eigenfunctions u^ix) = W123 (x) /W12 Ui(x) = 
Wi2i(x)/Wi2(x). A straightforward calculation shows: 

W M (x) = e -(T3+74)*( 73 _ lA )TW 12 (x + S)/W 12 (x) (19) 

where 

T 2 = (7? - 7 3 2 )(7 2 2 - 7 3 2 )(7! 2 " 7 4 2 )(7 2 2 " 7 4 2 ) (20) 

As one can easily see, W^(x) induces a second order Darboux displacement, V\{x) = 
Vq(x + 8), if and only if the following two numbers coincide: 



iartanh 

7i 

iartanh 

72 



7l(73 + 74) 



7i + 7374 

72(73 + 74) 



72 + 7374 



(21) 
(22) 



Notice that there are three regions in which (|19|) yields a non-singular Darboux transfor- 
mation ©, n x = {73,74 < 7l}> ^2 = {71 < 73,74 < 72}, ^3 = {73,74 > 72}, but only in 
f?2 can one achieve the consistency between ()21|) and ()22|). visualized by the intersection 
of two surfaces represented on Fig. 2. The points (73, 74) on the intersection curve provide 
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Figure 2: The intersection of two surfaces <|21|) and (|22|l provides the consistent data for the 
second order Darboux displacements of the 2-solitonic potential (|18|) . 

the precise data for the displacement, with 5 defined as the common value of (|21|) and 
(El- 

The approximate forms of all these phenomena can be observed for the "collage poten- 
tials" composed of fragments of the transparent wells. The simplest case is the periodic 
potential obtained by truncating 1)14(1 to a finite interval [—a, a] and then periodically re- 
peating all over R. One arrives at a periodic V{x) with T = 2a, for which the Schrodinger 
equation (^Q) can be explicitly solved for any E = k 2 , permitting to calculate the Lyapunov 
function 

l D = 1 t i^ri 1 ) ^ + - Mil cos2ka 

(with wo = 7otanh7oa). Notice that the equation defining the band edges \D\ = 2 is no 
more complicated than for the Kronig-Penney potential. The discrete energy level of the 
original potential (|14j) at Eq = —Jq belongs now to the lowest allowed band [HJ- The 
former ground level expands into the first spectral band and the interval (Eq,oo) splits 
into the subsequent gaps and bands. 







-1.5 




-a a 3a 



Figure 3: Second order Darboux transformations of the periodically continued one soliton po- 
tential with a = 5 and 70 = 0.9. (a) The initial potential; (b-c) The modified forms after using 
with iti, U2 chosen to be the pairs of Bloch functions for ot\ — —2 and 02 = —0.9. 

To illustrate the result of the Darboux transformation, we choose a = 5 and 70 = 0.9, 
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the numerical values of the lowest band edges becoming Eq = —0.8107, E\ = —0.8090, 
E v = 0.0001, E 2 = 0.1578, E 2 > = 0.1580; E 3 = 0.5926, E y = 0.5929. If now the 
transformation (jlj) is applied, the transformed potential V\(x) approximates very well a 
displaced copy of V(x): the effect looks as if the main body of Vq(x) was displaced, leaving 
only some tiny remnants in vicinities of the former peaks (Fig. 3). 



-0 .18 



-0 .33 

-a a 3a 

Figure 4: Second order Darboux transformations of the periodically continued one soliton po- 
tential with a = 2 and 70 = 0.4. (a) The initial potential; (b-c) The modified forms after using JIJ 
with ui, u-i chosen to be the pairs of Bloch functions for a% = — 10 and a 2 = —2. 

The phenomenon quite obviously imitates the behaviour of the original 1-soliton well, 
with an accuracy depending on the size of the repeated 1-soliton fragment (the bigger 
the fragment the better the effect). To examine this accuracy, we have performed the 
next computer experiment for the periodic potential composed of smaller pieces of (|14j) . 
taking 70 = 0.4 and a = 2. The lowest band edges are now Eq = —0.2664, E\ = 0.3204, 
Ey = 0.3817, E 2 = 2.1967, E 2 > = 2.2073, E 3 = 5.2838, E 3 > = 5.2885, and the results 
of the 2-order Darboux transformation are shown on Fig. 4. One can again observe a 
displacement affecting the main part of the original potential, though now the points of 
non-differentiability visibly resist the operation. 

An even more interesting effect occurs for the fragmented 2-soliton wells. By choosing 
the truncation borders ±a exactly at its minima, then repeating periodically and applying 
ijijl one gets a surprisingly exact picture of the displacement (Fig. 5). 
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-a a 3a 

Fi gure 5: Darboux operations on a periodically continued two soliton potential with 71 = 0.8, 
72 = 0.805 and a — 4.0279. (a) The initial potential; (b) the result of a 2-order Darboux trans- 
formation for a pair of Bloch functions with ou = —0.6, 02 = —0.0007 in the first energy gap 
(Ei = —0.6359, E v = 0). Notice a very good approximation to a finite displacement. 
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The effect once again cannot be perfect (the points of discontinuity of the third deriva- 
tive are now fixed), but the difference between the supersymmetrically transformed and 
displaced potential on our scale is practically invisible. 

Significantly, the existence of this kind of non-local phenomena (exact or approximate) 
tells a lot about the nature of the periodicity defects which can be supersymmetrically 
generated. To see this, let us consider a class of periodic potentials where the (non-local) 
displacements appear in their purest form. 



5 Supersymmetric transformations of the 
Lame potentials 

We refer to the Lame potentials, frequently considered in crystallography: 

Vo(x) = n(n + l)m sn 2 (x\m), n G N (23) 

where sn(x|m) is the standard Jacobi elliptic function. The potentials (j23[) have exactly 
In + 1 band edges, n + 1 allowed and n + 1 forbidden bands, and there exist analytic 
formulae for the eigenfunctions at the band edges (see e.g. j!7l). 

We shall show that the global SUSY displacement is a typical phenomenon in the 
subclass of Lame functions with n = 1, already for the first order Darboux transformations. 
In fact, we have 

Theorem 2. The Bloch eigenfunctions u^{x), u 1 /^(x) for a factorization energy 
a < Eq used in the first order Darboux transformations generate a displacement V\(x) = 
Vq{x + 8) of an arbitrary periodic Vq(x) if and only if 

U P{x)u l/f3 {x + 5) = c (24) 

where c is a constant. 

Proof. Suppose that u^(x) induces the Darboux displacement Vi(x) = Vq(x + 8). 
The standard Darboux theory and u^(x + T) = (5u^{x) imply that u l ^{x) oc l/u^(x) 
is the Bloch eigenfunction for hi with the same factorization energy a, so that u 1//3 (x + 
T) = (l//?)n 1///3 (x). As Vi(x) = Vq(x + 8), the coordinate change x — > x + 5 in the 
Schrodinger equation for u l /^{x) and the fact that the Bloch eigenfunctions are unique up 
to a multiplicative factor imply that u l ^{x) oc u l /^{x + 8). Hence 

uWlx + 8) = -4r^ => u l3 (x)v}IP(x + 8) = c (25) 
iiP{x) 

Conversely, suppose that the Bloch eigenfunctions of ho for the eigenvalue a satisfy 
u^(x)u 1 ^(x + 5) = c. Then perform a first order Darboux transformation using u^(x). In 
terms of u@(x) and ^^(x) the initial potential is given by: 

V (x) = 1 / +o= ; ; +a. 26 

Quite similarly, the final potential can be expressed in terms of u l l@(x) oc l/u@(x) = 
u 1 ^{x + 8)/c: 

Vi(x)= „ l/0( . +a= +a. (27) 
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By comparing and (f2*T|) one immediately sees: V\{x) = Vq(x + S). I 

Notice now that the criterion (I24[) indeed holds for the Lame potentials with n = 1. 
In order to prove that, consider the Bloch functions associated to the corresponding Lame 
equation (see Section 23.7): 

uf3 = ajxp + iJ) <j(x + a + u;') c _ c(a)(x _ XQ) 
a(xQ + a + oj') a{x + u') 

a{XQ — a + oj') a(x + oj') 

where xq is a fixed point in [0, T = 2K) selected so that n^(xo) = u 1 '^{xq) = 1, j3 = 
exp[2a£(w) — 2w£(a)], and u = K, u' = iK' are the real and imaginary half-periods of 
the Jacobi elliptic functions, a and £ are the non-elliptic Weierstrass functions, and the 
factorization energy a and a are related by 

a = |(m + l)-p(a), (29) 



where p is the well known Weierstrass function jzo]. Thus: 

U^(Xj« "'(Z + Oj — CT(xo+a+w / )CT(xo _ a+w /) CT (x+ w ') CT (a;+o+a;') e W 

By taking a = 5 we arrive at: 

u [x)u (x + d)- ^ + s + _ 6 + ujf) VI) 

where the right hand side does not depend on x, as was to be proved. It is not difficult to 
check that our criterion is not satisfied for the Lame functions with n > 1 ( |26| . Section 
23.7). 

Let us remind that the former results permitted to induce 5 = T/2 |71IH1E|- This is 
now recovered for a = Eq and for the unique Bloch eigenfunction u[x) = tpo(x) satisfying 
Mx)Mx±T/2) = c. 
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Figure 6: The translational effect of the 2-nd order Darboux transformation, (a) The initial 
Lame potential with n = 1 and m = 0.5; (b) The 2-susy equivalent. The factorization energies 
a.\ — 1.1, «2 = 1-4 belong to the first energy gap [E\,Ey), the displacement 8" = 0.747 ^ T/2. 
The final effect is very simple but it is not reducible to the nonsingular 1-st order steps. 



The existence of the corresponding second order displacements in the upper spectral 
gap follows from the inverse spectral theorems. Indeed, if the transformation (jlj) involves 
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the Bloch functions (|11I12() . it does not change the spectrum of Vq. As the shape of the 
Lame potential with n = 1 is uniquely defined by its band structure (cf. the inverse 
scattering theorems |321l33llM] . see also |22| p. 299, Theorem XIII. 91. b), we infer that the 
non-singular Darboux transformation of Proposition 1 can cause no more but a coordi- 
nate displacement. Once again, the abstract argument can be supported by the explicit 
formulae. In fact, by superposing a pair of Darboux operations © with the complex 
transformation functions Ui(x) given by (|28j) . where a = 5i are the complex displacements 
5i = 5 + u>', 82 = S'+u', one returns to the original Lame potential displaced by 5" = 5 + 5' 
(Fig. 6). The corresponding second order transformation has the real Wronskian: 



where 




T 2T 



Figure 7: A non-trivial result of the periodicity preserving second order Darboux transformation, 
(a) The original Lame potential with n = 3 and m = 0.5 (b) The Darboux deformed version. The 
factorization energies a\ = 2.15, «2 = 4.05 belong to the energy gap (Ei,Ei>), The global 
displacement affects the minima and maxima but simultaneously the potential is deformed. 

To examine the limitations of the method we have applied the second order Darboux 
transformations [To 1 l2"H l l%5 1 l3f) |l37| IH%] to the case n = 3 (see Fig. 7). A simple comparison 
with previous results shows that the effect is global but essentially new (one sees the 
displaced minima and maxima but also a non-trivial deformation). As it is already known, 
the first order Darboux transformations cannot displace the Lame potentials with n > 1 
(see |21|). The problem whether the same effect can be produced by higher order Darboux 
operations is still open (compare the discussion of Khare and Sukhatme 8 with Dunne 
and Feinberg [Jj). 

Once clarified the mechanism of the Darboux invariance, we have used the transforma- 
tion functions in (J3HU) defined as nontrivial linear combinations of the Bloch basis. Notice 
that the corresponding Darboux operations must affect the periodic structure, since on 
both extremes x — > ±00 the transformation function reduces to two different Bloch func- 
tions, causing two opposite asymptotic displacements. On Fig. 8 we show the periodicity 
defect and the injected localized state of the Lame potential (|23|) with n = 1, due to the 
1-st order transformation © with a < Eq. The Fig. 9, in turn, shows a defect of the same 
potential caused by the 2-nd order (irreducible) transformation Q), which has injected a 
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pair of localized states at the factorization constants a± = 1.2, «2 = 1-3. In both cases, 
the effect looks local but it is not (in fact, it is enough to compare minima and maxima 
of the initial and transformed potentials). The resulting periodicity defects arise as if a 
detail of the lattice was crushed by two opposite supersymmetric displacements, creating 
a "Darboux model" for the contact effects. 





(a) (b) 

Figure 8: The result of the first order Darboux transformation applied to the Lame potential 
with n = 1 and m — 0.5: (a) the supersymmetrically generated periodicity defect. Notice the 
asymptotic translational invariance; (b) the energy bound state for a — 0.35 injected into the 
infinite forbidden band (— oo,E ). 



To summarize, let us underline again an intriguing role of the translational invariance, 
which can appear either as an exact or approximate and/or asymptotic effect. This is an 
exceptional situation when the most elementary symmetry transformation, the unitary but 
nonlocal finite displacements, can be implemented by Darboux transformations typically 
producing local but non-unitary effects, a phenomenon specially convenient for the exact 
description of the contact effects (tentatively, including tiny quantum wells 39 ) . 






(a) (b) 

Figure 9: The result of the second order Darboux transformation (@J with the factorization 
constants a\ = f.2, ct2 = 1.3 applied to the Lame potential with n = 1 and m = 0.5: (a) the 
supersymmetrically generated periodicity defect. Notice the asymptotic translational invariance; 
(b) one of the energy bound states for a± — 1.2 injected into the forbidden band (E\, Ey). 
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